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Abstract. We study the baryon-quark phase transition in a molecular dynamics (MD) of quark degrees 
of freedom at finite temperature and density. The baryon state at low density and temperature, and the 
deconfined quark state at high density and temperature are reproduced. We investigate the equations of 
state of matters with different u-d-s compositions. Then we draw phase diagrams in the temperature- 
density plane by this simulation. It is found that the baryon-quark transition is sensitive to the quark 
width. 
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1 Introduction 

It is expected that quarks are observed as deconfined states 
at extreme environment: at high density and/or high tem- 
perature, where baryons or thermal pions overlap with 
each other, disappearance of hadron boundaries may give 
rise to the quark gluon plasma (QGP). Also the change of 
coupling constants and the quark mass due to the nature 
of the quantum chromodynamics (QCD), is thought to be 
one of the origins of the QGP [Q. Recently many theo- 
retical calculations are attempted to draw a QCD phase 
diagram 0EUE]. The critical temperature of the hadron- 
quark phase transition is predicted to be around 150 MeV 
by lattice simulations 0E]- For a finite density system, 
the bag model predicts in a qualitative study [7] that 
the critical density lies at a few times the normal nuclear 
saturation density. Many efforts are devoted to experi- 
mental search for the QGP by using high energy accel- 
erators. There, many indirect signals of the QGP have 
been observed but a definite conclusion has not 

been obtained since theoretical studies are not enough to 
characterize the properties of the QGP and hadron gas 
at present. It is not clear how hadron matter changes to 
quark matter and how the physical property changes at 
the hadron-quark transition. 

The Lattice simulation based on the first principle is 
the most reliable method for QCD. It has been applied 
to high-temperature and low-density systems under some 
approximations. At high density, however, the complex 
fermion determinant, or the sign problem, makes the sim- 
ulation practically impossible |llj . One of the approaches 
for treating finite density systems is a mean field theory 
which makes the most of its ability for uniform systems. 



There are many works with the mean field treatment, 
which are incorporated with the chiral restoration by the 
Nambu-Jona-Lasinio model ^2jj the bag model picture 
[T5| . the soliton model picture and so on. The dy- 
namics of hadron-quark transition, however, are not con- 
sidered in these cases. In these circumstances, other ap- 
proaches which are feasible for treating the finite density 
and the dynamics are being awaited. 

The molecular dynamics (MD) has been successful in 
treating many-body nucleon systems. To describe the struc- 
ture and the dynamics of many-body nucleon systems, 
several MD models were proposed |15II16II17II18| and ob- 
tained some remarkable results. Advantages of the MD 
simulation are that one needs very few assumptions a 
priori and that a single model can be applied to various 
problems. To investigate properties of hadron and quark 
matter and the dynamics of transition between these two 
phases in a unified manner, MD is a natural framework to 
be attempted. In this paper we apply a model similar to 
quantum molecular dynamics (QMD) ^] to many-body 
quark systems, within a framework of the non-relativistic 
quark model. Some pioneering works of applying MD to 
quark systems were presented where transitions caused by 
increasing density and temperature were studied in the 
view point of many-body dynamics: in 19,22] the Vlasov 
equation and the Vlasov+MD approach were used to ge t 
the equation of state of matter. In [2] quark MD was ap- 
plied to heavy-ion collisions where the creation of particle- 
antiparticle pairs was taken into account. In 22] . the color 
gauge symmetry was treated exactly and the meson ex- 
change potentials were introduced. In the above studies 
quarks were treated as classical particles and there was 
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ambiguity in the definition of the ground state of the sys- 
tem. In the present paper, we propose a model which has 
less ambiguity to describe the ground state and to define 
the temperature. With this model we study the mecha- 
nism of the baryon-quark transition and draw the phase 
diagram together with the equation of state (EOS) for a 
wide range of baryon density and temperature. 

This paper is organized as follows. In section [21 we ex- 
plain basics of the model in detail. Simulation results for 
the three kinds of matters: ud (matter containing the same 
number of u and d quarks; corresponds to symmetric nu- 
clear matter), udd {d quarks twice the number of u quarks; 
neutron matter) and uds (the same number of u, d and 
s quarks; A matter) by using two kinds of quark width 
are given in section In section 0] we apply the model to 
finite temperature systems and present phase diagrams. 
Summary is given in section [5] 



2 molecular dynamics for quark matter 

2.1 Wave functions and cooling equations 

We start with the total wave function by a direct product 
of n — 3A single particle quark Gaussian wave packets 
in coordinate and momentum spaces and state vectors \ 
with a fixed flavor, a color and a spin orientation, 
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where n = 3 A [A is the baryon number), L the fixed 
width of wave packets, and Ri and Pj are the center of the 
wave packet of i-th quark in coordinate and momentum 
spaces, respectively. Instead of antisymmetrization of the 
total wave function, the fermionic nature of the system is 
phenomenologically treated by introducing the Pauli po- 
tential which acts repulsively between quarks having the 
same flavor, color and spin orientation |23| . The equations 
of motion for R^ and are given by the Newtonian equa- 
tions, 
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When we search for the ground state (energy minimum 
configuration) of the system, we solve the equations of 
motion with friction terms, which we call "cooling equa- 
tions of motion" , 
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where \i r and /x p are negative frictional coefficients. The 
cooling is performed until the particles stop (R, = 0). In 
the ground state, the momenta P 4 have finite values be- 
cause of the momentum dependence of the Pauli potential. 

In order to simulate the infinite system (matter) by 
using a finite number of quarks, we use a cubic cell with 
26 mirroring cells under a periodic boundary condition. 
The cell size is chosen to be 6 fm throughout this paper. 



2.2 Effective Hamiltonian 

The effective Hamiltonian consists of three parts as, 

H = H o + Vpauli — Ispur, 

where Hq is the original Hamiltonian expressed as 
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We employ the quark-quark interactions as follows, 
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m where ry = |r» — rj\ and I means a light flavor u or d. 
The color dependent interaction V^oior consists of the lin- 
ear confining potential with an infrared cut-off at 3 fm 
and the one gluon exchange potential [22] with Gell-Mann 
matrices A a . To include the antisymmetric effect for the 
matrix elements of color space, we use effective values of 

(Xi\^^\Xj) eS = 4< X ,|AfA-| Xj > m The color force is 
approximately canceled between a colored quark and a 
white baryon made of quarks with three colors locating 
near to each other. V mcson consists of a, u> and p meson 
exchange potentials which act among light flavor quarks 
(u quarks and d quarks), and the <f> meson exchange po- 
tential which acts among s quarks, where u 3 is the third 
component of Pauli matrices. We modify the a exchange 
potential of Yukawa-type with a small non-linearity pa- 
rameter e. This corresponds to a density dependent po- 
tential or a method in the relativistic mean field theory 
(RMF) where higher order terms in a field are introduced 
to reproduce the saturation property for the symmetric 
nuclear matter. For simplicity, we have written Eq. (gj in 
the form of an operator. In practice, however, the power 
by 1 — e of the a exchange potential is performed after 
evaluating the expectation value. The parameters in the 
meson exchange potentials are adjusted to reproduce the 
baryon-baryon interactions as described in Sec. [21 C. 

Lack of the antisymmetrization is compensated by us- 
ing the Pauli potential, 
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where qo, po and C p are parameters determined in Sec.El 
C. Antisymmetrized wave functions are not used because 
it takes CPU time proportional to the forth power of a 
particle number 1 1 fill 1 T j . On the other hand, MD without 
antisymmetrization needs CPU times proportional to the 
second power of the particle number. At present, simu- 
lation with antisymmetrization is not practically possible 
for the system with several hundreds of particles. Further- 
more, the way to antisymmetrize wave functions with a 
periodic boundary condition has not been established yet. 
We have to subtract the spurious zero point energy of the 
center-of-mass motion of clusters |17H18| . 
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where Rij = |R S ; — Rj| and Mi is the "mass number" 
of fragment to which the wave packet i belongs, which is 
the sum of "friendship" with other particles. To apply 
the model to various densities, we introduce a density- 
normalized parameters of cluster separation in the friend- 
ship as 
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where po means the normal nuclear density 0.17 fm 



2.3 Choice of parameters 

2.3.1 Quark model parameter and friendship 

We use the constituent quark mass m u ,d = 300 MeV for 
light quarks and m s = 500 MeV for s quarks throughout 
this simulation. For the color potential, 900 MeV/fm is 
used for the string tension K and 1.25 for the QCD fine 
structure constant a s , which are typical in quark models 

The parameters in the friendship are chosen as clq = 
0.3 fm and wq = 0.5 fm, so that the sum of the friendship 
Mi « 3 when the system clusterizes as baryons and Mi ss 
1 when quarks do not make clusters. 

The width of quarks L is the most important parame- 
ter in this model since it is directly related to the density 
at which baryon-quark transition occurs. We use two dif- 
ferent widths by considering the masses of isolated nucleon 
N and lambda particle A. To minimize the energy (mass) 
of a nucleon, L becomes 0.43 fm. However, the nucleon 
mass is too large (about 2400 MeV) with this value of 
L. The reason of this overestimation is that the ground- 



state kinetic energy per quark in a nucleon, 



3h z 




Fig. 1. Energy per baryon of free Fermi gas. The lines show 
the theoretical exact values of ud, udd and uds matter. The 
marks are the kinetic energies calculated by MD with only the 
Pauli potential. 



has a large value of 351 MeV, while the kinetic energy of 
a quark in a nucleon is roughly estimated to be 50 - 80 
MeV from the uncertainty principle. Therefore, we employ 
"effective" widths L eS in evaluating Ei (Eq. (6)) and T spur 
(Eq. (10)). The ground-state kinetic energy per quark in 
a nucleon becomes 65 MeV if we use L eS = 1.0 fm. 

Our first choice is to use L = 0.43 fm, L° ff d = 1.0 fm 

and L° — 0.8 fm. We call this choice as "width parameter 
(I)" . This combination still gives a slightly large value of 
the nucleon mass, M/v ~ 1500 MeV, but gives a proper 
difference of masses between N and A. 

Once the effective width L for kinetic energy is em- 
ployed, the nucleon mass does not have a minimum regard- 
ing L. Then the second choice is made to give a proper 
nucleon mass value, Mn ~ 940 MeV by changing L to be 
0.33 fm (but with the same values of L cS as before). In 
this case, the difference of masses between N and A can- 
not be set to a desired value. We call this choice as "width 
parameter (II)". 
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2.3.2 Pauli potential 

For the Pauli potential, the parameters C p ,qo,po are de- 
termined by fitting the kinetic energy to the exact Fermi 
gas at zero temperature. We determine these values by 
solving the cooling equation ||3J) for ud matter, udd matter 
and uds matter where only the Pauli potential is consid- 
ered (23! ■ Figure 1 shows the classical and non-relativistic 
kinetic energy of the Fermi gas. The lines indicate the 
values of the exact Fermi gas energy for ud matter, udd 
matter and uds matter. The squares show the kinetic en- 
ergy of each matter, whose values are obtained by using 
following parameters: we adopt go = 1-6 fm, p = 120 
MeV and C p = 131 MeV for light quarks and C p = 79 
MeV for s quarks. Here three different colors and two dif- 
ferent spins are assumed for each matter. It is seen that 
the difference between ud matter and udd matter is small. 
By introducing heavy s quarks, the Fermi energy of uds 
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matter becomes lower than that of ud matter. Note that 
the Pauli potential gives a spurious potential energy to the 
system, which should be renormalized into other effective 
potential terms |2Sj- One possibility to avoid this prob- 
lem is, instead of using the Pauli potential, to maintain 
the Pauli principle of the system by stochastic rearrange- 
ments of particle momenta |25| . This model was originally 
developed for nucleon systems and was applied also to the 
quark system very recently |26j. 



2.3.3 Meson exchange potentials 

Though the meson exchange potentials do not affect the 
baryon-quark transition, we adjust the meson exchange 
potentials to discuss the EOS of strange and normal nu- 
clear matters. Prameters are determined to get appropri- 
ate ground state energy of "baryon matter" by frictional 
cooling with a constraint that quarks form baryon clusters 
(baryonization constraint). First we randomly distribute 
baryons which are composed of three quarks and solve the 
cooling equation with the baryonization constraint as, 
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here {i} means a set of three quarks in a baryon to which 
i-th quark belongs [22] • We call this cooling "baryon cool- 
ing" . In Fig. 2 baryon density dependence of the energy 
per baryon is shown for ud matter (corresponding to sym- 
metric nuclear matter), udd matter (neutron matter) and 
uds matter (A matter) . Note that in this calculation with 
the baryonization constraint, the color-dependent interac- 
tions are exactly canceled between the white baryons. For 
the a exchange potential, g aq — 3.09, m a — 400 MeV, 



1.2 fm and e = 0.1 are used and g u 
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Fig. 2. Binding energies per baryon for ud, udd and uds matter 
calculated by the baryon cooling. The saturation is seen for the 
ud matter, but it is not seen for the udd and uds matter. 



m w — 782 MeV and L u = 0.7 fm for the w ex change 
and g pq = 9.0, m p = 770 MeV and L p = 1.2 fm for 
the p exchange potential. In order to reproduce the well- 
known properties of matter, we have introduced the ef- 
fective widths for each meson exchange term. We fit the 
binding energy of ud matter to 16.5 MeV/nucleon at lpQ. 
The EOS of A matter, which depends on the 4> exchange 
potential is unsettled yet [22] ■ In order to obtain the sat- 
uration of uds matter, we need to introduce an attractive 
K exchange interaction between u-s or d-s quarks. How- 
ever the K meson exchange is prohibited according to the 
RMF model. Here we have only a repulsive (f> exchange 
interaction between s quarks for simplicity. The relevant 
parameters are g$ q = \/2g^ q = 7.04 |27|. = 1020 MeV 
and L,/, = 0.7 fm. Even if we use another value, e.g. ~ 1 
fm, the behavior of the EOS is not changed noticeably. 



3 Results for the finite density system 

Here we investigate the stability and the structure of quark 
matter in the ground stat e (zero temperature) for a wide 
range of density. We solve the cooling equations J2Jl with 
quark degrees of freedom. We call this cooling "quark cool- 
ing" in contrast to the baryon cooling. Snapshots of the 
ground state of uds matter with the width parameter (I) 
(i=0.43 fm) are displayed in Fig. 3. Three quarks in red, 
green and blue located within a distance d c iuster are con- 
sidered to be confined as a baryon 
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(15) 



where we use cluster = 0.7 fm. Quarks judged to be in a 
confined state are plotted in white and those deconfincd in 
its own color. At the normal baryon density lpo all quarks 
are confined as baryons. As density increases, some quarks 
become deconfined and the confined and the deconfined 
states coexist around 2.5po- Most of the quarks are decon- 
fincd at high baryon density p > 3po- The ground state 
energies of ud, udd and uds matters with the width pa- 
rameter (I) are shown by solid lines in Fig. 4. The vertical 
axis indicates the energy per baryon. For comparison, the 
case of baryon cooling is shown by thin dashed lines. The 
energy by quark cooling agrees with that by baryon cool- 
ing at low density. This means that quarks are confined 
as baryons. At a certain density, a baryon-quark transi- 
tion occurs and the energy decreases. The mechanism of 
baryon-quark transition in our dynamical model is as fol- 
lows. If confined quarks are released from a baryon, their 
kinetic energy, most of which has been the zero point en- 
ergy, decreases. On the other hand, the potential energy 
increases due to the confinement force. At above a certain 
density, the increase of potential energy at deconfinement 
gets smaller due to the existence of many other quarks in 
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Fig. 3. Snapshots of uds matter with L=0.43 fm at various 
densities. Using the criterion of Eq. p5^. confined quarks are 
plotted with white and deconfined quarks its own color. At lpo 
all quarks are confined. As the density increases deconfined 
quarks or its colors begin to appear. Most of the quarks are 
deconfined at 3po- 



the environment. In this way, the decrease of the kinetic 
energy is superior to the increase of potential energy. 

In our result with the width parameter (I), the critical 
density is lower than that usually expected. The ground 
state energies by the width parameter (II) (L=0.33 fm) 
for ud and udd matters is shown in Fig. 5. The critical 
point of the baryon-quark transition is higher than that 
given by the width parameter (I). This implies that the 
transition occurs due to the density overlap of a quark 
with those in other baryons. 

Figure 6 shows the energies for the three kinds of mat- 
ters under the electric-charge neutral condition with the 
width parameter (I). The difference compared to Fig. 4 
is that the energy of relativistic electron is added to the 



2000 



> 1800 
§ 1700 



§ 1600 




Fig. 4. Density dependence of energy per baryon for ud, udd 
and uds matter in case of width parameter (I) (L=0.43 fm). 
The dashed lines indicate the cases of baryon cooling and solid 
lines quark cooling. 
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Fig. 5. The same figure as in Fig. 4, but with the width pa- 
rameter (II) (L=0.33 fm). 




Fig. 6. Same as Fig. 4, but under the electric charge-neutral 
condition. 
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energy of ud matter. Below 6po the energy of udd mat- 
ter is the lowest. Judging from the behavior of the curves, 
however, ud matter becomes more stable than udd matter 
above 6pQ. 

The stability of strange matter has long been an at- 
tention of many nuclear and astrophysicists [281129] . Our 
present result shows that uds matter is not favored even at 
high density. However, we cannot give a definite statement 
since our EOS of uds matter has ambiguity caused by ig- 
norance of the A- A interaction. The width of s quarks and 
interactions concerning s quarks influence the stability as 
well. 

Here we have examined these three kinds of fixed flavor 
matter for simplicity. For more realistic studies of matter 
realized in nature such as the core region of neutron stars, 
however, a simulation of matter in beta equilibrium is nec- 
essary. The most stable matter with non-integer u-d-s pro- 
portion would be close to udd matter at lower density, and 
between udd and ud matter at higher density. With the 
width parameter (I) there are density regions where the 
slope of E/A versus p/po is negative in Fig. 6. Generally 
such a situation is not realized by forming non-uniform 
structure. In our case, however, the small size of the simu- 
lation box and the introduction of density-dependent pa- 
rameters lfF3*)l caused this problem. 

Now let us define a quantity Y which indicates the 
degree of deconfinement as, particle feels. 



Y = 



^ Kolor (Rj — R) 




1/2 



(16) 



where R is the position of a test particle in any color, 
V co \ OI is the color potential which the test particle feels and 
() R means the average over R which is sampled for 10000 
points. If all quarks are confined in compact baryons, the 
value of Y becomes zero. Figure 7 shows the density de- 
pendences of Y. The value of Y is small at low density. 
From 2po to 3po which corresponds to the density region of 
baryon-quark transition, Y increases gradually. This grad- 
ual change of Y does not necessarily indicate the second 
order phase transition, since the mixed phase in the first 
order phase transition may also show the similar behavior 
of physical quantities. 



4 extension to the finite temperature 

4.1 Introduction of an effective temperature 

We extend above MD to finite temperature systems. First 
of all, the definition of a temperature commonly used in 
MD is modified from the common expression due to the 
momentum-dependent Pauli potential according to Ref. 

EDI, 
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A popular method to control temperature in MD is known 
as Nose- Hoover method The equilibrium is achieved 
by introducing an imaginary heat bath which interacts 
with a real system through an artificial coordinate and 
its momentum. The extended Hamiltonian for thermostat 
systems is written as 



-Hnosc — H + 2<3£ 2 



3nT sot log it; , 



(18) 



with £ = P w /Q, where w is an artificial variable, P w is a 
momentum conjugate to w and Q corresponds to its mass. 
T set is a target temperature of the system which is given 
as a parameter. The equations of motion are 
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(20) 
(21) 
(22) 



Here H^ osc is conserved but the original effective Hamilto- 
nian H fluctuates. The value of T e g also fluctuates around 
Tset by the thermostat friction £. 
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Fig. 8. Averaged binding energy per nucleon for ud matter with L=0.43 fm at 0.3po (left), 0.5po (center) and lpo (right) are 
shown with solid line. The dashed line denotes the case in which quarks are constrained to form baryons. The arrows indicate 
the critical temperature. 




4.2 Results for the finite temperature system 

We investigate the thermal property of quark matter for 
a wide range of density and temperature. Figure 8 shows 
temperature dependence of the averaged energy per baryon 
for ud matter at 0.3po: 0.5po and Ipo with L = 0.43 fm 
from the left by solid lines. Simulation results with the 
baryonization constraint are also shown by dashed lines for 
comparison. A baryon-quark transition can be seen as the 
change of specific heat. We define the critical temperatures 
by the maximum-specific-heat points and are indicated by 
arrows in this figure. The critical temperature becomes 
small as density increases. The lines without the bary- 
onization constraint (normal calculation) disagree with 
those with the baryonization constraint because quarks 
can move in the baryon in case of normal calculations. 
The same quantity but for uds matter is shown in Fig. 9. 
The critical temperatures for uds matter are 5-15 MeV 
larger than that for ud matter at the same density. This is 
due to the presence of the s quarks having a larger mass 
than u and d quarks. 

We draw a phase diagram by using our results of ud 
matter for L=0.43 fm in Fig. 10 and for L=0.33 fm in Fig. 
11. Baryon phase is realized at lower side of density and 
temperature, and quark phase is realized at higher side in 
both diagrams, but baryon phase in Fig. 11 is larger than 
that of Fig. 10. This means that the quark width affects 
the phase transition due to the increase of temperature as 
well as the increase of den sity. 



5 Summary 

Quark many-body systems were studied by MD where the 
ground state was defined in a definite manner in terms 
of a Pauli potential. The EOS were reproduced for three 
kinds of baryon matters with effective meson exchange 
interactions between quarks. Baryon-quark transition is 
seen when baryon density increases. We have used two 
quark widths, L = 0.43 fm and L — 0.33 fm. The density 
at which baryon-quark phase transition occurs is differ- 
ent for the two widths since the transition is caused by 
the overlap of the quarks. In case of the larger value of 
L = 0.43 fm the baryon-quark transition occurs at rather 
low density, around 3po. For L = 0.33 fm case, the tran- 
sition occurs around 5pg, which is consistent with other 
theoretical calculations. 

We have compared three kinds of matters with differ- 
ent u-d-s compositions: ud, udd and uds. Our results show 
that ud matter is the most stable among them and strange 
uds matter is not stable. For the EOS of uds or A matter, 
there would be still room for improvement to take into 
account the future experiment. For the color interaction, 
the color magnetic interaction is necessary as a color and 
flavor dependent interaction. The medium effects of con- 
stituent quark masses by the chiral symmetry restoration 
are also important to discuss the stability. A possibility to 
include the medium effect is to use the density dependent 
quark mass model (DDQM) (22]. However, the quark mass 
is derived there from the global density and is common for 
all quarks. This prescription may be too simple since the 
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Fig. 10. Phase diagram with the density and temperature for 
ud matter in case of L—0.43 fm. 
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Fig. 11. Phase diagram with the density and temperature for 
ud matter in case of L=0.33 fm. 

individual particle motion is essential in MD. A similar 
discussion can be done for the coupling constant a s . 

We have extended our model for the finite temperature 
systems. The baryon-quark transition can be seen as a 
change of specific heat. The critical temperature for uds 
matter is higher than ud matter by 5 - 15 MeV due to 
the heavy mass of s quark. Like the critical density, the 
critical temperature gets higher in case of L=0.33 fm than 
that of L=0.43 fm. 

It is necessary to include qq creation/annihilation pro- 
cess for high temperature. Lack of antisymmetrization and 
the asymptotic freedom are also open problems. 

Y. A. is grateful to T. Tanigawa, H. Koura, T. Tatsumi, T. Hat- 
suda, Y. Maezawa and T. Endo for their valuable discussions. 
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